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Experimental evidence for Majorana bound states largely relies on measurements of the tunnel¬ 
ing conductance. While the conductance into a Majorana state is in principle quantized to 2e 2 //i, 
observation of this quantization has been elusive, presumably due to temperature broadening in the 
normal-metal lead. Here, we propose to use a superconducting lead instead, whose gap strongly 
suppresses thermal excitations. For a wide range of tunneling strengths and temperatures, a Ma¬ 
jorana state is then signaled by symmetric conductance peaks at eV = ±A of a universal height 
G = (4 — 7t)2 e 2 /h. For a superconducting scanning tunneling microscope tip, Majorana states 
appear as spatial conductance plateaus while the conductance varies with the local wavefunction 
for trivial Andreev bound states. We discuss effects of nonresonant (bulk) Andreev reflections and 
quasiparticle poisoning. 


Introduction .—Motivated by possible applications in 
quantum information processing [1, 2], topological super¬ 
conductors hosting Majorana bound states are currently 
under intense investigation [3-5]. Based on the super¬ 
conducting proximity effect, various realistic platforms 
have been proposed to support Majorana states includ¬ 
ing topological insulators [6, 7], semiconductor nanowires 
[8, 9], and atomic chains [10-16]. Although these systems 
are available in the laboratory, the experimental observa¬ 
tion of unique Majorana signatures remains challenging. 

A widely employed diagnostic tool is the tunneling con¬ 
ductance of normal metal-superconductor junctions, in 
which Majoranas manifest themselves as characteristic 
zero-bias peaks [17, 18]. Experimental signatures con¬ 
sistent with theoretical predictions have been observed 
in quantum wires [20-22] and atomic chains [23, 24], 
However, it is a major challenge in these experiments 
to uniquely distinguish Majoranas from conventional 
fermionic subgap states. Spin-polarized subgap states 
such as Shiba states bound to magnetic impurities [25-28] 
or Andreev bound states in a magnetic field can exhibit 
a zero-energy crossing as a function of exchange inter¬ 
action or Zeeman energy [29-31]. Thus, such fermionic 
states may accidentally occur at zero energy and give 
rise to similar conductance features. As magnetic im¬ 
purities or external magnetic fields are also required for 
the most relevant realizations of topological superconduc¬ 
tors, such trivial conductance peaks can generally not be 
disregarded. 

In contrast to fermionic subgap states, Majoranas ex¬ 
hibit a celebrated quantized zero-bias conductance of 
2e 2 //i [17-19]. Unfortunately, this has so far proved dif¬ 
ficult to observe in experiment. The Fermi distribution 
in the metal lead is smooth on the scale of the temper¬ 
ature T, which strongly limits the experimental energy 
resolution. When temperature is larger than the tunnel 
coupling, the Majorana peak is broadened and the zero- 
bias conductance is reduced. Even at low temperatures 
(e.g., T = 60 mK in Ref. [20]), it may be difficult to ob¬ 
serve the quantized peak height as multichannel effects 


limit the relevant tunneling strength [32]. Quasiparticle 
poisoning may also lead to deviations from quantization. 
A fermion-parity breaking rate exceeding the tunnel cou¬ 
pling broadens the peak and reduces its height. This 
requires one to work at temperatures below the lowest 
fermionic excitations in the topological superconductor. 

In this paper, we show how a robust conductance sig¬ 
nature of Majorana bound states can be obtained by em¬ 
ploying superconducting leads. In striking contrast to 
normal-state contacts, effects of thermal broadening are 
strongly suppressed for a superconducting lead because 
quasiparticle excitations are exponentially suppressed ~ 
exp(—A/T) by its superconducting gap A. Majorana 
bound states no longer appear as zero-bias anomalies but 
rather as two symmetric peaks in the differential conduc¬ 
tance G = dl/dV which occur when the BCS singularity 
of the superconducting gap lines up with the Majorana 
bound state, i.e., at the thresholds eV = ±A. These 
peaks have a universal height 

G m = (4-7t) 2 ® 2 , (1) 

which persists over a wide range of tunnel couplings. 

This yields particularly striking evidence when employ¬ 
ing a scanning tunneling microscope (STM) with a super¬ 
conducting tip which allows for spatially resolved mea¬ 
surements. This has previously been used to map out 
bound state wavefunctions in conventional and unconven¬ 
tional superconductors [23, 24, 33-36]. Here we propose 
that such maps can clearly distinguish between Majo¬ 
ranas and trivial zero-energy bound states. Indeed, the 
peak conductance is uniform in the vicinity of Majorana 
states and a conductance map exhibits a characteristic 
mesa or plateau structure. In contrast, the conductance 
of trivial subgap states exhibits a spatial pattern which 
is governed by the bound-state wavefunction. 

In addition, STM measurements allow for systematic 
studies as a function of tunneling strength by varying 
the tip height. It was recently demonstrated [37] that 
this can be exploited to probe quasiparticle relaxation 
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processes. In the present context, varying the tunnel¬ 
ing strength may help to identify Majorana signatures 
despite competing effects such as nonresonant Andreev 
reflections or quasiparticle poisoning. 

Subgap conductance for Majorana bound state.- At 
subgap voltages eV < A + A s and zero temperature, the 
tunneling current between superconducting tip or lead 
and substrate (with gap A s ) flows by multiple Andreev 
reflections. Near the threshold e\V\ = A, the differential 
conductance dl /dV is dominated by single Andreev re¬ 
flections from the sample. For tip locations far from the 
zero-energy bound state in the sample, this yields the 
familiar peak in dl/dV due to the singular densities of 
states of incoming electrons and outgoing holes. In the 
vicinity of the bound state, tunneling is further enhanced 
by the zero-energy resonance [37-39]. 

Formally, the subgap current due to single Andreev 
reflections from the sample can be expressed as [40-42] 

I{V) = 4e7r 2 t 4 J ^Tr[G eh {r,u)Gl h (r,u)} 
x p(u-)p(w+)[n F (uj-) - n F (u + )], (2) 


in terms of the electron and hole tunneling rates 
r e //j(w) = 27rf 2 |((| 2 p(w= F ). While the integrand in Eq. (4) 
has a resonance denominator, its behavior is peculiar due 
to the strong energy dependence of the tunneling rates. 
Specifically, the square-root singularity of the BCS den¬ 
sity of states implies that the integrand involves a charac¬ 
teristic energy scale uj t = (M 2 v 0 \((r)\ 2 A/2) 2 / 3 which 
depends on a fractional power of the tunneling rate from 
a normal tip y n = 27rf 2 i'o|£(r)| 2 . In the weak-tunneling 
regime < A, we can write 

iMJ-f Jt -A- T , (5 ) 

for 0 < rj <C A, where g = eV — A measures the voltage 
from the threshold A. In the vicinity of the threshold, 
i) < W(, the resonance denominator is dominated by the 
second term and we obtain I(V ) = (4 — n)(2e/h)(eV — 
A)9(eV — A) and thus Eq. (1). The entire peak lineshape 


d I_ 
dV 


(4-tO^A 


eV- A 


( 6 ) 


where t is the amplitude for tip-substrate tunneling, 
u>± = uj ± eV, uf(uj) denotes the Fermi function, and 
the superconducting tip enters through its BCS density 
of states p(oj) = vo0(\u\ — A)|w|/Vw 2 — A 2 with the 
normal density of states at the Fermi energy. Spin or 
subband degrees of freedom are accounted for by a pos¬ 
sible matrix structure of the anomalous retarded Green 
function G e h(r,Lo) of the substrate at the tip position r. 
In terms of its Lehmann representation, G e h(r,co) has 
contributions from both the bound state and the above¬ 
gap continuum. In the following, we first consider the 
resonantly enhanced Andreev current from a Majorana 
bound state and subsequently discuss the contribution of 
the quasiparticle continuum. 

For e|Vj ~ A, we can approximate (<*;_)—rip (w+) ~ 
sgnF in Eq. (2), up to corrections of order exp(— A/T). 
This insensitivity to temperature is a key advantage of 
superconducting leads. The bound-state contribution to 
the substrate Green function is 


G(r,w) 


WHV’I r) 
ui + iT/2 


( 3 ) 


Here, (r\if) = [£(r), ±©C(c)] T denotes the local 

Bogoliubov-de Gennes wavefunction of the Majorana 
bound state with 0 the time-reversal operator. The 
broadening T = 2i (ip\Y,\ip) of the bound state is induced 
by the tunnel coupling to the lead. The corresponding 
self energy E = —*7rt 2 diag[p(w_),p(w + )] is diagonal as 
Andreev reflections in the lead can be neglected near 
e\V\ = A. 

Inserting Eq. (3) into (2) yields (for V > 0) [37, 43] 


r -l 


dui 


r e (w)r ft ,(w) 


w 2 + [r e (w) + Thito )} 2 /d 


( 4 ) 


involves the function A(x) which vanishes for x < 0, 
jumps to A(0 + ) = 1, and falls off with a small negative 
differential conductance tail at large x, cp. Fig. 1. 

Thus, the differential conductance between a conven¬ 
tional superconductor and a Majorana state exhibits a 
peak which is independent of tunneling strength and Ma¬ 
jorana wavefunction. While the peak height is close to 
the quantized Majorana peak height 2e 2 /h for a normal- 
metal lead, there are several differences: (a) There are 
two symmetric, finite-bias Majorana peaks at eV = ±A 
rather than a single zero-bias peak, (b) the conductance 
peak is strongly asymmetric with a discontinuous step at 
the threshold, and (c) the width of the peak is set by ui t 
with its sublinear dependence on junction transparency. 

The threshold discontinuity in the conductance per¬ 
sists even when including the contributions of the quasi¬ 
particle continuum in the substrate Green function. To 
see this, we model the substrate superconductor by a 
2x2 Nambu Green function g(u>,r). For a topolog¬ 
ical substrate, this is appropriate for perfect spin po¬ 
larization (spinless p -wave superconductor). Including 
the tunnel coupling to the tip through the self energy 
E as given above, the substrate Green function becomes 
G = g[ 1 — Eg] -1 . We first focus on the vicinity of the 
bound state where the conductance is dominated by An¬ 
dreev reflections from the bound state. By straightfor¬ 
ward calculation and expansion of g in u> [40], we find 


G(r,u) 


(r\if){if\r) 
u> — A(w) + /2 


( 7 ) 


This differs from the pure bound-state contribution by 
the additional term A(w) = 7r 2 f 4 wdetg(w,r)p(w_)/9(w+) 
in the denominator which involves the determinant (in 
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FIG. 1. (color online) (a) Differential conductance vs bias 
voltage near the threshold eV = A for Majorana (solid line) 
and Andreev state with |u| = |n| (dashed line). For a Majo¬ 
rana, the conductance exhibits a step of height (4 — n)2e 2 /h 
at the threshold. For an Andreev state, the conductance has 
a smooth onset, cf. Eq. (10). Both peaks have a negative- 
differential conductance dip at high voltages. Inset: Graph 
of f(x) as defined in the main text, (b) Spatial conduc¬ 
tance maps for Majorana (left) and Andreev state (right) for 
Wt(0)/JA = 5. The Majorana gives rise to a conductance 
plateau whereas the Andreev state exhibits a pattern reflect¬ 
ing the spatial dependence of the ratio ( u/v ) 2 . The Majorana 
conductance drops far from the bound state when the broad¬ 
ening exceeds u>t- 


particle-hole space) of the bare substrate Green func¬ 
tion. While the determinant of the bound-state contri¬ 
bution to the Green function vanishes, this is no longer 
the case when including the quasiparticle continuum. At 
subgap energies away from bound states, the Green func¬ 
tion g(co, r) is a hermitian 2x2 matrix, so that det g(uj, r) 
and hence A (ui) are real. Thus, we find 


3/2 n ' 

“V ) 2 

/rj+uJ 

( 8 ) 

For a Majorana state, the real part of the resonance 
denominator must vanish exactly at u = 0. Indeed, 
particle-hole symmetry further constrains det g(ui,r) to 
be an even function of w which can be approximated 
as a constant at small oj (see [40], where this conclu¬ 
sion is confirmed by model calculations). Then, we find 
A(w) oc t^cj/y/r] 2 — ui 2 near the threshold. Even with this 
term, the denominator in Eq. (8) remains dominated by 
the divergent tunnel broadenings ~ aS’J' 2 /\Jr] ± u and the 


‘M 


(F) = x/ 


du> 


-v vV-w 2 ( w _ A )2 + (- 


3/2 


discontinuous conductance step as well as the universal 
value of the threshold conductance in Eq. (1) persist. 

In experiment, the square-root singularity of the BCS 
density of states of the tip may be broadened intrinsi¬ 
cally due to higher-order processes or effectively due to 
experimental resolution. The universal threshold con¬ 
ductance persists as long as u h exceeds this broadening. 
This condition also determines the spatial extent of the 
conductance plateau, r < 4£ln[w t (0)/<5A]/3, where £ is 
the Majorana localization length, ujt( 0) denotes the value 
of tot. at the center of the Majorana bound state, and S A 
is the broadening of the tip density of state, cf. Fig. 1(b). 
Of course, a well-resolved Majorana peak also requires 
cut <C A s , i.e., the tunnel broadening needs to be small 
compared to the induced gap. If the peak is not fully 
resolved, it is suppressed below the universal value and 
its height may vary as a function of space. 

For tip locations far from the bound state, the tunnel¬ 
ing conductance is dominated by conventional (‘“nonres¬ 
onant”) Andreev reflections. These still yield a threshold 
peak due to the singular tip density of states in F e and 
Th, but are not enhanced by a bound-state resonance. 
For a one-dimensional p- wave superconductor, this con¬ 
ductance peak has height ~ 1.3 Gm and width ~ AT 2 
quadratic in the junction transparency T oc t 2 [40]. Ob¬ 
serving the conventional Andreev peak thus requires that 
the broadening of the tip density of states is small com¬ 
pared to ~ A T 2 . This is a much more stringent condition 
than for the resonant Andreev peak as the width of the 
bound-state peak u> t oc t 4 / 3 involves a lower power of t. 
We note that in a typical STM experiment [37], conven¬ 
tional Andreev peaks can be resolved only for small tip- 
sample distances, while bound-state signatures persist to 
much weaker tunnel couplings. 

Subgap conductance for Andreev bound state. —These 
results should be contrasted with those for trivial zero- 
energy Andreev bound states. For concreteness, con¬ 
sider an s-wave superconductor with conserved spin [44] , 
whose Bogoliubov-de Gennes description decomposes 
into two independent spin sectors that interchange un¬ 
der particle-hole transformations. A zero-energy An¬ 
dreev state corresponds to two Bogoliubov-de Gennes 
wavefunctions, {r\i/j+} = [it(r), w(r)] T and (r\i/j-) = 
[0u(r), — 0u(r)] T , one in each sector. An analogous cal¬ 
culation [40] yields the threshold current 

I A (V) = 2I M (V)f(\u(r)\ 2 /\v(r)\ 2 ). (9) 

Reflecting the two zero-energy wavefunctions, the max¬ 
imal threshold conductance is twice that in the Majo¬ 
rana case, Ga = ZGm, and realized for the particle- 
hole symmetric case |ri| = |v|. In general, the peak 
conductance depends on the ratio of electron and hole 
wavefunction at the tip position. This dependence 
is captured by the dimensionless function f{x) = 
dzy/l — z 2 j{x\J\ — z + \J\ + z ) 2 which takes on 
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values between 0 and 1 and is plotted in Fig. 1(a). The 
function satisfies f{x) = /(1/a;) as the two spin sectors 
contribute equally. In the limit of large particle-hole 
asymmetry, Ga ~ GMmin(|u/i;| 2 , \v/u\ 2 ) <C Gm • The 
lineshape of the conductance peak is similar to the Ma- 
jorana peak, with a width of order io t upon replacing £(r) 
by max{ u(r), v(r)}. 

Our results imply that the height of the conductance 
peak allows for a clear distinction between a conventional 
Andreev bound state and a Majorana state. Even when 
f(u 2 /v 2 ) ~ 1/2 for one location of the STM tip, mov¬ 
ing the tip to another location modifies the conductance 
peak height for a conventional bound state, tracking the 
ratio of electron and hole wavefunctions. In contrast, the 
conductance map exhibits a characteristic mesa struc¬ 
ture for a Majorana state, see Fig. 1(b). In non-STM 
tunneling experiments, changes of parameters (e.g., gate 
voltages) which affect the Majorana wavefunction should 
leave the peak height unchanged for a Majorana but not 
for a conventional Andreev bound state. 

As there is no locking of the bound state to zero en¬ 
ergy, also the continuum contribution is distinctly differ¬ 
ent for conventional Andreev states. The two spin sectors 
are described by separate 2x2 Nambu Green functions 
which map into one another under particle-hole transfor¬ 
mations. This is quite unlike the Majorana Green func¬ 
tion which maps onto itself. For each sector, det g(io,r) is 
therefore no longer an even function of io and will gener¬ 
ally have a singular contribution oc l/u> at the threshold 
so that A(w) ~ T 2 A s A/^f — lo 2 . These general ar¬ 
guments can be confirmed explicitly for Shiba states in 
s-wave superconductors [40]. Near the threshold, the res¬ 
onance denominator in the expression for the current is 
now dominated by A(w). As illustrated in Fig. 1 by a 
numerical evaluation of the current, this suppresses the 
conductance step. Analytically, we find that just above 
the threshold, the conductance increases linearly, 

2p 2 1 pV — A 

Ga(V) ~ - A), (10) 

and matches with the conductance obtained from Eq. (9) 
for eV — A T 2 A. We note that this suppression of 
the conductance step depends on T and can thus be 
probed by varying the tip-sample distance in an STM 
experiment. This may serve as an additional signature 
to distinguish between Majorana and conventional An¬ 
dreev bound states. 

Effects of quasiparticle poisoning. —So far, we only in¬ 
cluded bound-state broadening by the tunneling contact. 
At finite temperatures, the bound-state occupation also 
changes by inelastic transitions to other subgap states or 
the quasiparticle continuum in the sample [45]. We ac¬ 
count for these processes by an additional contribution 
*r qp /2 to the self energy of the bound-state Green func¬ 
tion Eq. (7). This does not affect the Andreev current 
at the threshold, where the denominator is dominated 



FIG. 2. Total threshold conductance for a Majorana state 
(tot) along with the single-particle contribution (sp) as a func¬ 
tion of LOt- The single-particle contribution affects the con¬ 
ductance only a in a window of transmission values, where 
ujt ~ r qp . While the maximum is of order 0.2 Gm, the po¬ 
sition of the maximum in tunneling strength depends sensi¬ 
tively on temperature (through r qp ). Inset: Line shape of 
the total conductance as a function of voltage away from the 
threshold, for different ratios of c Jt/T qp . 


by the diverging tunnel coupling. However, the over¬ 
all weight of the peak is reduced by a narrowing of the 
linewidth by a factor (w t /F qp ) 2 once r qp > w t , see Fig. 
2 (inset). 

In addition, quasiparticle poisoning generates a single¬ 
electron current I s which involves tunneling of single 
particles followed by inelastic transitions from the zero- 
energy bound state to other bound states or the quasipar¬ 
ticle continuum [37]. For a Majorana state, we find near 
the threshold eV = A (with analogous results applying 
for Andreev bound states) [40] 


ts __ _f_ f j _ r qp [r e (o;) -i- rfa(cc)] _ 

m_ 4 hj u 2 + [r qp + r e (w) + r /l (w)] 2 /4' 

For weak and strong tunneling, this yields [40] 

qs ^ f (wt/r qp ) 3/ ujt <c r qp , 

^ 1 T qp /u>t LOt Tqp. 


( 11 ) 


( 12 ) 


Figure 2 shows that this single-particle contribution as¬ 
sumes a maximum of ~ 0.2Gm when ui t ~ r qp . However, 
it can be easily made negligible by tuning the system 
away from this maximum through varying temperature 
or tunneling strength. 

Conclusions. —We show that conductance measure¬ 
ments with superconducting leads constitute a promising 
technique to identify Majorana states. The presence of 
Majoranas is signaled by conductance peaks of universal 
height which are largely unaffected by thermal broaden¬ 
ing, a key obstacle in previous experiments with normal- 
metal contacts. We discuss strategies to systematically 
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rule out parasitic effects such as quasiparticle poisoning 
or trivial subgap states. The proposed setup is read¬ 
ily available in the laboratory and, in fact, has already 
been realized in previous experiments [23, 24, 46, 47]. 
(Notice, however, that temperature was comparable to 
the induced gap in the STM experiments performed to 
date, precluding observation of the universal conduc¬ 
tance, and that the nanowire experiments focused on 
zero-bias peaks.) Our results also imply that quasipar¬ 
ticle poisoning rates can be extracted from systematic 
measurements as a function of tip height and tempera¬ 
ture. 
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SUPPLEMENTARY MATERIAL 

I. GENERAL FORMULA FOR ANDREEV CURRENT INTO SUBGAP STATES 


In this section, we outline the derivation of the tunneling current in Eq. (2) of the main text. This standard 
calculation is included to make the presentation self contained and closely follows the derivation for 2 x 2 Green 
functions presented in Ref. [37]. We describe the tunneling contact by the Hamiltonian H = Hr + Hr + Hr, where 
the three terms refer to the lead (tip), the sample, and the tunnel coupling. The superconducting tip with chemical 
potential /i and gap A is described by the BCS Hamiltonian 


H l 


dk 

(27t) 3 


E& -Am cl ,k a + (a 4, ,kt^L,-14 



(13) 


where £& = k 2 /2m — /i and Cl, k<x (c' L k ) annihilates (creates) an electron in the tip with momentum k and spin er. 
The sample Hamiltonian generally takes the form 


H r = 


dx E $R, a { x YH-R,a<j> (x)cr^(x), 

a,a' 


(14) 


where 'H.R^a’ix) is the Hamiltonian in first quantization and cr iCT (x) annihilates an electron with spin a at position 
x in the sample. We choose the superconducting order parameters in tip and sample to be real such that the 
superconducting phase difference 4>{t) enters the tunneling Hamiltonian 


Ht{t) = X! [ iel0(T)/2 CL, CT (O,T)c^ CT (a;, 


-) + te ^ (t)/2 c], (T (x,t)c L!(T (0,t) 


(15) 


where r the time argument, t is the hopping strength, and Cl i<t (0,t) = f dkci,k,x(T)/(27r) 3 annihilates an electron in 
the tip at the tunneling contact, which is located at the origin. The sample is contacted at position x and we suppress 
the position arguments in the following for simplicity. The time-dependent phase difference between the tip and the 
sample, </>(r) = cj) 0 + 2 eVr, depends on the voltage V applied to the junction. 

We evaluate the current from the Heisenberg equation of motion / = — eN r = ie[Ni,, Hr), where Nr is the 
electron-number operator of the tip. Taking the expectation value, we obtain 

I ( T ) = | Tr { r z [K t )Gr L {t,t) - Gr R (t, T)t*(r)] } , (16) 

where r 2 is a Pauli matrix acting in Nambu space, i(r) = te lT ^ T ^ 2 T z , and we have introduced the lesser Green 
function in Nambu and spin space with matrix elements given by (G^)^ = i(\E , j 3j \E' Q ,j). Here the ’P a j are components 

of the Nambu spinor = ^c a> |, c ai j,, j_, —with a, /3 = L, R. 

We are interested in tunneling processes to lowest order in the tunneling amplitude and therefore neglect Andreev 
reflections in the tip, which give rise to multiple Andreev reflections. As these processes involve several single-particle 
tunneling events they enter only at higher orders in the tunneling amplitude. This approximation is exact when the 
sample is spin polarized (e.g., a proximity-coupled semiconductor nanowire in a strong magnetic field) and one of 
the two spin components is fully normal reflected. In this case, spin-flipping Andreev reflections in the tip do not 
contribute to transport. 

We follow the nonequilibrium Green function approach described in Ref. [37, 41] setting the offdiagonal elements 
of the tip Green function in Nambu space to zero and denoting the diagonal elements by ^l(w). We obtain the dc 
current 

1 = / duJ ^{ G> R ee ^)9U“ - eV) ^ G<’ e » 5 >(u; - eV) - g>(co + eV)G<’ hh (u>) + ff£(w + eV)G > R ’ hh (w)} , 

(17) 

where Gr(uj) = f dTidr 2 exp[iw(r 1 — t 2 )\Gr(ti, t 2 ) is the sample Green function in presence of the tip, (e, ft.) are 
indices in Nambu space denoting particle and hole components, and the trace is taken in spin space. This expression 
has been derived in Ref. [37] for the case when tip and sample are both spin-conserving s-wave superconductors. 
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According to the Langreth rule, the lesser Green function of the sample can be written as 

G%=g%+ SrErGr + 9R^G a R + g R T, R G R 

= (1 - g r R^ r R)~ l 9%( 1 + %) + G r T, r G r . (18) 

Similar expressions also exist for the greater Green function. The first term involving g R gives rise to a single-particle 
current which at subgap energies requires inelastic processes in the sample. The second term leads to the current 
carried by Andreev reflections. The sample self energy due to the presence of the tip can be written as 

£(w) = t 2 diag{g L (u]-),g L {uj-),g L (uj + ),g L {uj + )), = u ± eV. (19) 


We now focus on the Andreev current and defer the discussion of the single-particle current to Sec. V. Thus 
neglecting the first term in Eq. (18) we find Eq. (2) of the main text, 

I(V) = ~^J du} Tr {G T R h {u)g^ (w + eV)G R he (u)gf;(u} — eV) — G r R eh (cj)g^ (uj + eV)G R he (cv)g R (w — eU) j 

+ j du Tr + eV> >G r R e {^)gl{^ - eV)G a / h {u) - g<(uj + eE)G^ /te (w)^(w - eU)G^ e,i (w)} 

J duj \\G e R {u)\\ 2 p{u> + eV)p(uj — eV) [rif(u — eV) — n R (ui + eV)] (20) 


47r 2 et 4 


where G r R h is the electron-hole block of the (retarded) Green function of the sample and ||G|| = ^Tr (GGl) denotes 
the Frobenius norm of matrix G. Here we have used the relations g^igj) = 2Trinp(ui)p(uj) and g^{uj) = — 27ri(l — 
tif(u}))p(u>), where uf{oj) is the Fermi distribution function, p(u>) = |w| 6(u) 2 — A 2 )/-\/w 2 — A 2 and vq is the normal 

density of state at the Fermi energy in the tip. Due to the step functions in p(ui ± eV) the integration interval is 
restricted to u £ (—(eU — A), eV — A). Note that in this interval the self energy is purely imaginary. 


II. CONDUCTANCE FOR ZERO-ENERGY BOUND STATES 

In this section, we calculate the conductance for isolated Majorana or Andreev states at zero energy, cf. Eqs. (6) 
and (9) of the main text, neglecting the contributions of all other states in the sample. While the main text focuses 
on Andreev states in s-wave superconductors, we also consider more general spin structures here. 


A. Topological superconductor with Majorana bound states 


We first consider a topological superconductor substrate with a single zero-energy Majorana state to the tip. The 
Majorana wavefunction has the form d>o(x) = u l( x )*i u l( x ), ^ u t( x )*) T which maps onto itself under a particle- 

hole transformation. Neglecting contributions from other states, we can approximate the sample Green function by 
3 m(W) x,x) = <S>o(z)$ > o( x )/ UJ - Note that we suppress position arguments throughout this section. Including the 
coupling to the tip the full Green function can be written as 

G = (21) 

W — ZjM 

where Em = 4>jE$o is the self energy projected onto the Majorana bound state. We obtain 


G r eh H 


1 

w + int 2 |C| 2 [p(w+) + p(u>-)} 


u t u i 

uf 


u t ^ 

J 


( 22 ) 


where we introduced |C| 2 = |tt-|-| 2 + |wj,| . Thus we find 


l|G eh .H|| 2 


_jc£_ 

W 2 + 7T 2 t 4 |C| 4 [p(w+) + p(w-)] 2 


(23) 









and using Eq. (20) we obtain the current 


rfnrfH 
w 2 + (rfM + rfM)74 


(24) 


with = 27rf2 ICI 2 p 7 t ). 

We now evaluate the current near the threshold eV = A + 77 , 77 <C A. The current at the opposite threshold 
eV = —A follows from I(—V) = —I(V). At low temperatures, T <C A, we can set np(ui — eV) — np(uj + eV) ~ 1 and 
the Majorana current reads 


9 p r> 

Im(V) = y I <iu 

To lowest order in 77 , one can approximate 


rf( w )rf 7 ) 


o ^ + (rfM + rfM) 2 /4‘ 


This yields 


p(u> ± eV) ~ 


2 e 


rK 9(r) ± uj) 

2 y/rf^uj 


Im{V) zz “ 7r ) A ( ? ?M) 
h 


as given in Eq. ( 6 ) of the main text, where we have defined 


(25) 


(26) 


(27) 


A(x) = 


4 — 7T 


dz 


Z*X 3 + 


f // - + 1 ) 

Vv / l + 2 ) 


At the threshold we find 


A(0) = 


4 — 7T 


dz 


l + VT^l 


= l 


(28) 


(29) 


which yields Eq. (1) of the main text. At large voltages, 77 w t , we instead find a negative differential conductance 
dl/dV oc — 1 / 77 3 in agreement with the lineshape shown in Fig. (1) of the main text. 


B. Non-topological Andreev states at zero-energy 

A non-topological zero-energy Andreev bound state is characterized by two Nambu spinors 

*+ = (")■ *- = ( 9 e V „) < 30 > 

where u = (iq-, rt|) T , v = (tq, tq-) T are functions of space and 0 = —ia y K, is the time-reversal operator with K, the 
complex conjugation. The Lehmann representation of the real space Green function is thus a 4 x 4 Matrix in Nambu 
and spin space 


4> + <E>l_ + 4>_4>I 
5M = — 


u> 


= -( $+,$- ) 

UJ 


(31) 


When spin is a good quantum number the Green function may be reduced to a 2 x 2 Matrix in particle-hole space 
only. In the case of pure s-wave pairing, where Cooper pairs are formed from electrons with opposite spin, we can 
set ui = Vf = 0. The spinors $ + and 4>_ then belong to the orthogonal subspaces spanned by (c-j-, 0, ct, 0) and 

(0, cj,, 0, —cl). The Green function decomposes into two 2x2 blocks, which are related by particle-hole symmetry and 
have equal contributions to the current. In the opposite case of a spin polarized p -wave superconductor we can set 
u l(x) = v\.{x) = 0. Now 4> + and <h_ belong to the same subspace spanned by (c-f,0,0, — cj) and the Green function 
is a single 2 x 2 matrix. 
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It is therefore useful to first discuss a general 2x2 Nambu Green function 

_ (9ee 9eh 

9 ~ \9he ghh 


(32) 


We will return to the 4x4 case in Sec. IIB 3 when discussing superconductors with spin-orbit coupling. The coupling 
to the tip can be included through the self-energy S r (w) = — * 7 rf 2 diag(p(w_), p(uq_)). The Green function of the 
coupled system is obtained from the Dyson equation G = g{ 1 — Eg) -1 . We find 


G r eh M = 


g e h(w)w 


w - A(w) + incut 2 (g ee (u])p(uj -) + g hh {u>)p(u + )) 


where A(w) = ujn 2 t 4 p(co + )p(u}-) det g(cu). Using Eq. (20), we obtain the current 


m = 


An 2 et A 


r eV-A 


dui 


W 2 \g e h(u )\ 2 p{uj-)p{u>+) 


/A-eV 


[w - A(w )] 2 + n 2 cu 2 t 4 l(g ee (uj)p(uj-) + g hh (u])p(u>+)} 


2 ’ 


(33) 


(34) 


We now calculate the conductance for a zero-energy Andreev state in the limiting cases of pure s-wave and spinless 
p-wave pairing and for a general spin-structure of the order parameter in the presence of spin-orbit coupling. The 
result for the three cases are compared in Fig. 3. 


1. s-wave pairing 

When ui = v-f = 0 the Green function in Eq. (31) is block diagonal. The 2x2 block in the basis (c-j-, ct) reads 


9 = 


u t v i 




(35) 


We find detg(w) = 0 and thus A(w) vanishes. Using Eq. (34), we arrive at the current as given in Eq. (9) of the main 
text, 


2e r eV ~ A 

U{V) = -r / du 

n JA-eV 


t^mt^m 


w 2 + (r^M + r£M)74 


(36) 


where r^(w) = 2nt 2 |rt|| 2 p(w_), r^(w) = 2nt 2 p(u>+). We have included an extra factor of two to account for 
the second 2x2 block of the Green function which yields an equal contribution as a consequence of particle-hole 
symmetry. Near the threshold when eV = A + 77 , rj <C A we find 


lA(V) = 2I M (V)f(\u(r)\ 2 /\v(r)\ 2 ), 


(37) 


which can be obtained by a similar analysis as for the Majorana bound state in the previous section. The dimensionless 
function 


f(x) = - - / dz\J\ — z 2 /{xy/1 — 2 + V1 + z) 2 

4-7r 


takes on values between 0 and 1. Thus, the threshold differential conductance is 


dlA 

dV 


eV=A 


= ^f(\u(r)\y\v(r)\*)(4-n). 


(38) 


(39) 


2. p-wave pairing 


In a spin-polarized p- wave superconductor Cooper pairs we can set iq = iq = 0. The Green function of an Andreev 
state reduces to a single 2 x 2 matrix 


1 ( Kl 2 + M 2 2w Wt ^ 

lj y 2u^vf |u-f | 2 + |w-f | 2 y 


( 40 ) 
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We have introduced A(w) = 7 r 2 t 4 p(w + )p(w_)(|u-f \ z — \v-f\ 2 ) 2 /uj and using Eq. (34) we obtain 


Ia(V) = 


I 67 r 2 ef 


2 0 f 4 rv 


du> 


|u t u t | p(w)p(w+) 


UJ ~ 


p(u-)p(uj + ) 7 r 2 \t 4 \ 


(41) 


2 - H 2 ) + 7 r 2 t 4 (|u t | 2 + K| 2 ) 2 [p(w + ) + p(w _)] 2 


with eV = A + p. Close to the threshold 77 —>• 0 the denominator is dominated by A(w) oc 1 / 7 y 2 and we generically find 
/ oc rj^9(p). Hence the conductance G oc p 3 9(p) is continuous in contrast to the case of s-wave pairing. At certain 
points in space it may be possible that |u^| = |u-f| in which case the conductance still jumps at the threshold. Slightly 
moving away from such points should restore the smooth onset of the conductance at the threshold. The conductance 
exhibits a peak at the characteristic scale p ~ u> t with up = (max{|ii^.| 2 , |t'-y| 2 }^ 0 t 2 ^A) 2 ^ 3 as shown in Fig. 3. 


3. Generic case with spin-orbit coupling 


Realistic proposals of topological superconductors typically involve a mixture of s-wave and p-wave pairing. In 
particular, such pairing arises in any superconductor with spin-orbit coupling. In this case all components of u and 
v are generically nonzero. The full sample Green function including the coupling to the tip due to the self-energy in 
Eq. (19) can be written in terms of Dyson series 



w 


<b + ) 


*1 


E( 4>+ 



A straightforward calculation reveals 


G r _ 

eh ~ 


(w — Eh_)u + 57+(Cv) 


/t - 


(w- E; + )Cv + E;_u (Cu)t 


_ w (e; + + S r _) + e; + sl_ - e;_el + 


where the projected self-energies are 


E++(w) = iirt 2 [||u|| 2 p(w_) + ||v|| 2 p(w+)] 
E!1_M = iirt 2 [||v|| 2 p(w_) + ||u|| 2 p(w + )] 
S+_(w) = -wrf 2 (u, 0 v) [p(w_) + p(w+)] 
£_+M = -I?rf 2 ( 0 v, u) [p(w_) + p(w+)] = -E+_(w)*, 


(42) 


(43) 


(44a) 

(44b) 

(44c) 

(44d) 


and (•, •) is the inner product. In the above derivations, we have used the anti-unitarity of the time-reversal operator 
and that 0 2 = —1, namely (u, ©v) = (v, 0^u) = — (v, @u). Then the norm can be written as 


\\G eh \\ 2 = W 


{w 2 


2u > 2 + Wn 2 t 4 Y {[p(u >+) 2 + p(cj-) 2 ]Z + 4 p(oj + )p(uj-)Y} 

IE 7 T 2 t 4 [(p(w + ) 2 + p(co-) 2 )Y + Zp(uj + )p(u}-)]} 2 + IEw 2 7 T 2 f 4 (Z + 2)[p(u>_) + p(uj + )} 2 


(45) 


where 


W = u 2 v | 2 + (u, 0 v) 2 

(46a) 

= (l|u|| 2 ||v || 2 - |(u, 0 v)| 2 ) /w 

(46b) 

(ll u || 4 + |v | 4 — 2 (u, 0v) 2 ) /W. 

(46c) 


The differential conductance is then a function of Y, Z and p/up, where uif = Av^ii 2 t 4 W. The parameter Y 
interpolates between s-wave pairing (Y = 1, Sec. IIB 1), where the threshold conductance is maximal, p-wave pairing 
(Y = 0, Sec. IIB 2), where the threshold conductance is zero. In Fig. 3 we show the conductance for different values 
of this parameter. 



















11 



(eV-A )/d)f 

FIG. 3. Conductance near the threshold. The parameters for the s-wave pairing case are Y = 1, Z = 4.25. For the p-wave 
pairing case, Y = 0, Z = 1.13. And in the case with spin-orbit coupling, Y = 0.05, Z = 0.74. 

III. CONTINUUM CONTRIBUTION TO THE CONDUCTANCE 

In this section, we provide details on evaluating the contribution of the quasiparticle continuum to the conductance 
near the threshold. We identify two different effects of the continuum states, namely, (a) a possible smoothening of 
the step in the bound-state conductance [see Fig. 1(a) and Eq. (10) of the main text] due to interference between 
resonant and nonresonant Andreev reflections and (b) additonal contributions to the conductance from nonresonant 
Andreev reflections. The latter effect becomes important when probing the conductance away from the bound state, 
where resonant and nonresonant Andreev reflections compete as discussed in the main text. 


A. Continuum effect on resonant Andreev reflections 

1. General considerations 


As follows from Sec. II, the lineshape of a zero-energy Andreev state in an s-wave superconductor closely resembles 
that of a Majorana state. In particular, the conductance exhibits a step at the threshold eV = A in both cases. We 
now revisit these two cases and discuss whether this step is robust when the quasiparticle continuum is taken into 
account. As discussed in the main text, we find a suppression of the conductance for the Andreev state close to the 
threshold, while the step remains robust for the Majorana. This behavior is reminiscent of tunneling from a normal 
metal tip, where Majorana states appear as a robust zero-bias conductance peaks while zero-energy Andreev states 
generically exhibit zero conductivity. 

For simplicity, we focus our analysis on samples in which spin is a good quantum number and that can be described 
in 2 x 2 Nambu space. An example of a more general model is discussed in Sec. IV. In particular, we consider an 
s-wave superconductor with a zero-energy Andreev state and a spinless p-wave superconductor with a Majorana state. 
In these cases the Lehmann representation of the Green function reads 




$o(aO$o(z 


E 


dE 


$E,n( X )®E,n( X ) 


|B|>A a 


UJ — E 


(47) 


where the spinor <h 0 describes a single zero-energy bound state and e, u are continuum states above the gap A s with 
energy E and index n labeling the degeneracy. Using this expression in the Green function in the presence of a tunnel 
coupling in Eq. (33) and expanding to lowest order in energy, we obtain 

$o(r)$o(r) 

w — A(w) + i(r e + Y h )/2 


G r eh M 


(48) 
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To determine the effect of the extra term A(w) oc det < 7 ( 0 ;) we analyze the energy dependence of det < 7 ( 0 ;) in more 
detail. The function g(uj) has a simple pole at to = 0 and branch cuts on the real axis for |w| > A s but is analytic 
elsewhere. The residue of g(uj) at w = 0 is which has zero determinant and we can thus write 

lim uj 2 det g(uj) = 0 (49) 

Col—>-0 

and expand det g(ui) as a Laurent series around w = 0 

OO 

detg(w)= ^ c n uj n . (50) 

n =—1 

Furthermore g(uj) is hermitian at nonzero subgap energies and thus its determinant is real. We now determine the 
lowest-order cotribution to det g for Majorana and Andreev states. 


2. Majorana states 


The Green function of a topological superconductor satisfies particle-hole symmetry Ucg(uj)U' c = —g*(—uj) with a 
unitary operator Uc■ This yields det g(uj) = — det g*(— uj) and reality further requires detg(w) to be an even function 
of uj. Hence, we obtain c_i = 0 and A (uj) ~ uit 4 p(oj + )p(oj-) 7 which approaches a constant in the limit ij = eV —A — > 0 + 
with |w| < g. The Majorana contribution to the current reads 


Im(V) 


2 e r rfHCH 

h J Q (UJ - A ) 2 + (Tf (w) + Tf H) 2 /4' 


(51) 


The rates TW diverge at the threshold and the continuum term A(w) becomes negligible. The threshold conductance 
of Majorana states thus remains unaffected by the continuum states. 


3. Zero-energy Shiba state 

In contrast to Majoranas, the 2x2 Green function describing zero-energy Andreev states generically does not 
satisfy particle-hole symmetry and thus c_i ^ 0. As an example of a system with trivial zero-energy Andreev state, 
we consider a magnetic impurity in an s-wave superconductor, which induces a bound state localized at the impurity, 
known as a Shiba state. The Hamiltonian describing a Shiba state localized at the origin due to a magnetic impurity 
can be written in first quantization as 


Hs(x) = Hbcs(x) + ( Vt z - JSa z )5(x), 


(52) 


where and eq are Pauli matrices in Nambu and spin space, where V and JS are the potential scattering and 
exchange coupling strength. Since Hs is block diagonal in spin space, we only need to deal with one of the blocks, say 
cr z = 1. The other block related to the first by particle-hole symmetry contributes a second (equivalent) zero-energy 
state as mentioned above. The sample Green function gg at the impurity position has the form [37] 


SsM = 


2 uja — (1 — a 2 + /3 2 )-y/A 2 — uj 2 


uj + (a + /3) V / A 2 - uj 2 A s 

A s ui + (a - /3) V / A 2 - w 2 


(53) 


where v s and A s are the normal density of states and gap of the sample and we introduced the dimensionless 
parameters a = ttvqJS and /3 = ttuqV. The Shiba bound state energy is given by 


to — A s 


l~a 2 +/3 2 


(54) 


yj(l — a 2 + /3 2 ) 2 + 4a 2 

Thus, since we focus on zero-energy bound states we set a 2 = 1 + /3 2 . After a straightforward calculation, we obtain 


wdetg(w) ~ 


2 a 


(55) 
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for small uj. From the Green function we also obtain (for V > 0) 

(56) 

(57) 


M 2 = 7ns s A s (a + Jo? - 1), 
|u| 2 = F v s A s (ol — 'Jo? — 1). 


Near the threshold, for small 77 = eV — A, the term ex A(w) in the denominator of Eq. (34) dominates and we find 


Ia(V) 


2 e 16 ? 7 2 a 2 
X z/ 2 !/27r 3 f 4 A' 


(58) 


Thus the conductance (11a/ dr) is zero at eV = A and then rises linearly as described by Eq. (10) of the main text. 
This originates from the interference between resonant and nonresonant Andreev reflection, which suppresses the 
differential conductance at eV = A. In order to estimate the width of this suppression we determine the value of 77 
when the term A in the denominator of Eq. (34) becomes of the same order of magnitude as the tunneling broadening 
~ ujt 2 [g ee p(u>-) + g hh p(ui+)], i.e., 


i/q At 4 w det g(ui) 
V 



which yields 


V 


vlv^A 



(59) 


(60) 


B. Nonresonant Andreev reflection 

We now determine the conductance due to nonresonant Andreev reflections from the quasiparticle continuum in 
the substrate. To evaluate the conductance contributions from the continuum, we focus on a specific model of a Dirac 
Hamiltonian with a domain wall, where the mass changes sign. We first calculate the real space Green function and 
then evaluate the conductance as a function of separation from the impurity. 


1. Topological superconductor Green function 

As an example of a topological superconductor hosting a Majorana state we calculate the Green function of a Dirac 
Hamiltonian with a domain wall 


HDirac = -iv F d x T x - m[9{x) - 9(-x)]t~ (61) 

in the spinless Nambu basis (ip,if'). Here v F is the velocity of the Dirac fermion and m is the effective mass. The 
system satisfies particle-hole symmetry {C,HE>; rac } = 0 with the charge conjugation operator C = t x JC, where /C is 
the complex conjugation. This is a generic model for the low-energy behavior of a topological superconductor close to 
the phase transition. Specifically, a spinless p-wave superconductor with a chemical potential close to the bottom of 
the band can be approximated by Eq. (61). To compute the eigenstates of this Hamiltonian it is convenient to solve 
for the rotated Hamiltonian e * WT 'y/ 4 7 p Dirace - l7rT !// 4 . Notice that the charge conjugation operator is now given by t z IC. 
One readily shows that the domain wall supports a single Majorana state at zero energy with the wavefunction 

*-<*> = (0 • < 62 > 

All other states in the system are extended with energies above the gap to. The wavefunctions of the continuum 
states can be determined using the ansatz 

iM*) = ( T } e) ) eiP " m ~ X) + m{x)] + ( Tv S FP±e) ) e ~ iPX {C6{ ~ X) + D9{x)] ’ (63) 
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where v F p = V£ 2 — rri 2 and e is the positive eigenvalue. The matching condition at x = 0 leads to 
A + B + C + D = 0, v F p{A - B - C + D) =f e(A - B + C - D) = 0. 

For positive-energy solutions, we get 

A + D 


A + C vpp y 

For incoming waves from x = —oo we require D = 0 and obtain 


/mV 

\VF PJ 


_ V FP € A^ n 


_ _ VFP 


Similarly, for outgoing waves to x = —oo we have A = 0 and thus 

This yields the scattering states 


QOUt VFP pout 


j^out _ _ 6 T VfP jjout 


(64) 

(65) 

( 66 ) 

(67) 


€ n (x) = A 

r e ut (x) = D 


—m9(—x) — m^r-O(x) 


(e - vfp)O(-x) - H££( e _ VfP )0(x) 
e + vfP I -to 
vfP — e 


vfP — e 
e 


—TO 

vfP + e 


0(—x)e 


9{x)e ipx + 


m9(x) — m^ E 9(—x) 

( vfP + e)9(x) + v F p— f ±£ 6*(—a;) 


0 —ipx 


( 68 ) 


One can easily check that incoming and outgoing states are orthogonal {4>° ut \ l l ) T') = 0. From the normalization 
conditions \\ipl n \\ — IIV’e“ t || = 1 we obtain the coefficients A = [2Le(e — v F p )] -1 / 2 and D = [2Le(e + v F p )] -1 ^ 2 for a 
system of size L. The Green function can now be obtained from its spectral decomposition 

>f(a;)^“(x , ) t , {C^{x)){C^{x')f' 


g(x,x';w) = 


tpoixWoix'Y 


CJ 


E 

e>0, 

a=in,out 


OJ — € 


lo + e 


(69) 


where we have used the charge conjugation operator C to relate the negative solution to the positive ones. Focusing 
on x = x' = 0 for simplicity we find the wavefunctions of the continuum states 


VC(o) = - 


l 


v F p 


a/2 Le(e - v F p) e 
Using Eq. (69), we obtain the Green function 


m 

e — vpp 


4C\ o) = 


1 


v F p 


ff(°j°; w ) — 2vf4v V f 

Here we have used 


:£ 


e>0 




i = —— ( 

(o 

ij 

2 vf ix \ 


a/2 Le(e + v F p) 


—i 


—TO 

e + VfP 


— (w/to ) 2 


jh _ I dP — 

T.2^1 ~ / Ott _ 


L 


e>0 


27r 2ttvf 


de 


•/e 2 — to 2 


(70) 


(71) 


(72) 


Notice that the matrix commutes with r y , and thus commutes with the rotation operation introduced at the beginning. 
We find the Green function for the original Hamiltonian i?Dirac 


m ft \ rn / y/1 - ( uj/m ) 2 

5(0, 0;w) = 7rzy s — I 

w \ —* 


a/1 - (w/to) 2 


(73) 


where = (27 tuf) _ 1 is the normal density of states. 

Away from the mass domain wall, the sample Green function at x > 0 can be computed similarly from Eq. (68) 
and (69), which yields 


5ee(*r) = ghh(x) = 


o-2z/« _ 


V 1 - ( co/m ) 2 V W TO/ 

i - e - 2:c /£ 


9eh(x) = ghe(x)* =TTV s [i — e 2x/i - - . • 

V W \/l - (w/to.) 2 / 


(74) 


where £ = [27 t^ s TOa/1 — (w/to) 2 ] 1 = vf/tu^/I — (w/to) 2 is the coherence length that characterizes the exponential 
decay of the Majorana wavefunction. 
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FIG. 4. Differential conductance of the topological superconductor with a mass domain wall. The three curves correspond 
different tip locations measured from the domain wall. For all curves we have set AT 2 /m = 0.016. 


2. Differential conductance away from the Majorana state 


The Green function has a pole at zero energy whose amplitude decays exponentially away from the domain wall 
reflecting the probability density of the Majorana bound state. This singular part fully determines the threshold 
conductance and yields the quantized value Gm- Notice that lim^^o w det g = 0 guarantees the robustness of this 
quantization consistent with the general argument based on particle-hole symmetry given in Sec. Ill A 2. The con¬ 
ductance at voltages above the threshold depends on the competition between the singular part and a nonsingular 
contribution which originates from the quasiparticle continuum. The continuum contribution can be readily evaluated 
for x —> oo, where the Green function describes a homogeneous p -wave superconductor, 


9 M 


\J TO 2 — iff 1 


uj m 
m uj 


(75) 


Using this expression in the current formula in Eq. (34) we find 


/ = 


47r 2 e |t| J 


~ T a7 " 2 


duj 


dz- 


2 2 2 
p+p-ir v^m 


(m 2 — Ul 2 ) (1 + 7T 4 I/ 2 t 4 /9 + p_) 2 + 7T 4 l/ 2 t 4 W 2 (/9 + + p_) 2 




> 2 > 


-i (2Vl - z 2 +AT 2 /vY 

ticti 

(277/ATVT^ 


(76) 


where we introduced the transmission probability of the junction T = 7 r 2 ^ s i/ot 2 and assumed AT 2 <C m. This yields 
the conductance 

r*l 


G(V) = 


8 G 


M 


dz 


4 - tt 7 0 “ (l + (2r ? /Ar 2 )vT^ 2 ) 3 ' 


(77) 


The conductance has a maximum at 77 ~ AT 2 with a magnitude ~ 1.3 Gm and a peak width ~ AT 2 . In Fig. 4 
we show the differential conductance vs voltage for different separations between tip and domain wall. As the tip 
is moved away from the domain wall, nonresonant Andreev reflections become more important and the conductance 
peak becomes narrower and slightly higher. 


IV. DOMAIN WALL AT THE QSHI EDGE 

While the quantized conductance step for the Majorana state derived in the main text is model independent, its 
robustness to continuum effects was only shown for spinless models described by a 2 x 2 Hamiltonian. As an example 
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of a more realistic model, we consider the edge of a quantum spin Hall insulator (QSHI). We calculate the threshold 
conductance including the continuum contribution for the case when the sample is contacted at the domain wall. We 
find the universal step at the threshold to be robust in line with the results for the spinless case. This suggests that 
the conductance step at the threshold may be robust even for more general models of topological superconductors. 

The edge of a QSHI, with a domain wall at the origin, can be described by the first quantized Hamiltonian 

H(x) = - iv F d x a z T z + A (x)r x + B{x)a x , 

A(*> = M - 


where 


m(x) = m [9{x) — 8(—x)]. 


Applying the unitary transformation U = e * e 

l i T y, we have 

( M 0 

0 

-iv F d x \ 

u^u = unu ] = 

0 

—m( x) 

-iv F d x 

m(x) 

0 

0 

—iv. F d x 

0 


\ -iv F d x 0 

0 

-M / 


(78) 


Let M m, and focus on the inner block, which decribes a Dirac field with a domain wall created by a mass jump, 
as discussed in Sec. IIIB 1. According to Eq. (73), the Green function at the domain wall is 


<?(0, ui) = 7 TU S 


UJ 


(o o 
0 
0 


°\ 


\0 


—l 

0 


1-S o 


o 7 


(79) 


We rotate it back to the original basis g(0,w) = Ug(0,oj)W and obtain the full Green function 

G(0,w) = g(0,w)(l - £(u%(0,w)) _1 , 


(80) 


where the self energy is given by 


£(w) = t 2 dia.g(g{u-),g(u-),g(u+),g(uj + )). 


(81) 


Calculation reveals 
G r ’ eh { 0,w) = 


'Kv s m 


w [1 + i 4 7T 4 i/|(p(w_) + p(u+)) 2 } + 2iTT 2 u s mt 2 -Jl — ^ [p(w_) + p{u+)} 

* yj 1 - £ ^ in (^f) t 2 ( P (u>-) + p(w+)) \ 

^l-^-in(^)t 2 (p(uj-)+p(uj + )) -i ) 


Denoting T e / h = 4-7r 2 z v s mt 2 p(uJ T ) we find 


G eh { 0,w)f = 


2(7 TVs ) 2 [2w 2 - w 2 (l - (T e + T ft ) 2 /16m 2 )] 
w 2 [1 - (T e + T h ) 2 /16m 2 ] 2 + (T e + T h ) 2 /A 


(82) 


We readily verify that the denominator is dominated by the second term at the threshold. In particular, we see 
the continuum contribution in the denominator appears proportional to w 2 as for the spinless case. Hence, the 
conductance quantization at eV = A is unaffected by the continuum contribution in this model. 
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V. SINGLE-PARTICLE CURRENT CONTRIBUTION TO THE CONDUCTANCE 

In this section we analyze the single-particle tunneling current that can flow in addition to the Andreev current 
discussed above. A single-particle current is possible if the quasiparticle occupying the bound state can relax to 
the quasiparticle continuum. Such relaxation can occur, e.g., due to inelastic transitions assisted by phonons or 
photons. We concentrate on the Majorana state, although a similar analysis is possible for a trivial Andreev state. 
We neglect the quasiparticle continuum and include relaxation as a phenomenological parameter, without specifying 
its microscopic origin. We first derive the single-particle current [given in Eq. (11) of the main text], from which we 
then calculate the threshold conductance from single-particle tunneling [Eq. (12)]. 


A. Self energy due to relaxation processes 


Let us assume that the local environment of the Majorana state introduces relaxation processes, induced by phonons 
or photons, via the self-energy E p h. The substrate Green function can be determined from the Dyson series 

9R = 9 + g^phg + <?Sp h <?Ep h <? + ■ • •, (83) 

where g is the Green function without relaxation processes. We approximate the bare substrate Green function g by 
the Majorana contribution, 

g(co) = \iPm)-(4>m\, (84) 

UJ 

and project the self-energy onto the Majorana subspace introducing 

r qp = 2Im(^ M |Eph(0)|V’M), (85) 


where we approximate the self-energy by its value at zero energy. The invariance of the Majorana state under particle- 
hole transformation guarantees that the expectation value ('0jvf|E p h(O)|'!/ , M) is purely imaginary. Thus, the retarded 
and advanced Green functions of the Majorana state read 


r.a ( \ 

gA M 


IV'mXV’mI 

uj ± iTqp/2 


( 86 ) 


In quasi-equilibriunr, we can express the greater and lesser Green function in terms of the retarded and advanced 
Green functions, 


= /M(s&( w ) - g r R (u)) = 2 + 

OJ T 1 qp/^ 

E > (0) 

9r{u) = -(! ^ /H)(5fi(w) - s£M) = ^2^2 / 4 \iPm){iI>m\, 


(87) 

( 88 ) 


where f(u) is the quasi-equilibrium distribution function and we used the relations 


- *E< h (0) = T qp /, *Ep h (0) = T qp (1 - /). (89) 

These terms are the rates for emptying and filling of the delocalized fermion formed by the Majorana at the contact 
and a second one far away. Note that local transitions can change the occupation of this state. Since this fermion has 
zero energy the two rates are equal j'E^O) = — iE^O) = r qp /2 according to detailed balance. 


B. Expressions for the single-particle current 


We can now evaluate the current in Eq. (17). Besides the Andreev current in Eq. (20), we find the single-particle 
current from the first term in Eq. (18) 


Im(V) = ^Jd“~ 


r qp [T e (u)n F (uj-) - r h (u})n F {u>+)] - r qp [r e (u;)(l - n F {uj-)) - T h (a;)(l - n F (ui+))] 


w 2 + (T e (w) + T^(w) + Tqp) /4 
which gives Eq. (11) in the main text if we take n F (ui-) ~ 1 and n F (u> + ) ~ 0 assuming T <C A. 


( 90 ) 






18 


C. Analysis of current 


We compute the single-particle current near the threshold at eV = A + 77 with r\ <C A 


Im(V) ^ 


2 h 


a+ v w t 3 / 2 r qp ( 

t 4P \ y/rj—uj y/rj+uj ) 

duj ---— 


— A—rj 


OJ* 


Uf /2 e(ri+Lu) u>t /2 8(ri—u) 


/V+u 


/r)-ui 


„ 3 /2-p 

euj t r qp 


duj 


/UJ 


( W - d) 2 + ^ + 2 


3/2 


2 > 


to lowest order in 77 . We thus obtain the differential conductance 

r , M _ 2e 2r qp f°° du} 
u m( x ) — 


huit 


{u-x) 2 + (u 1/2 +r qp /(2w t ))' 


I 2 ’ 


eV- A 


(91) 


(92) 


Now we focus on the conductance at the threshold, namely x = 0. For weak tip-substrate tunneling, u> t <C r qp , we 
find 


G 


S 

M ~ 


2e 2 r 

i QP 

huj t 



__ 

w 2 + (r qp /( 2 w t )) 2 


27re 2 wf^ 2 
h p 3 / 2 ’ 

1 qp 


(93) 


where the 2 ;-integration is elementary. In the opposite limit of strong tip-substrate tunneling, tot r qp , we can 
neglect the contribution of r qp in the denominator. In this limit, we find the peak conductance 


G s m ( A) 


r diU — 

h^t Jo [cu 2 + 1 / (d] 2 

27re 2 r qp 
9 h u>t 


( 94 ) 


The parametric dependence of G S M is summarized in Eq. (12) of the main text. 



